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■ Abstract 

We construct a novel approach, based on thermodynamic geometry, to 
characterize first-order phase transitions from a microscopic perspective, 
through the scalar curvature in the equilibrium thermodynamic state space. 
Our method resolves key theoretical issues in macroscopic thermodynamic 
constructs, and furthermore characterizes the Widom line through the 
^ ■ maxima of the correlation length, which is captured by the thermody- 

namic scalar curvature. As an illustration of our method, we use it in 
conjunction with the mean field Van der Waals equation of state to pre- 
dict the coexistence curve and the Widom line. Where closely applicable, 
it provides excellent agreement with experimental data. The universality 
of our method is indicated by direct calculations from the NIST database. 
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Macroscopic properties of matter undergo discontinuous changes along a first- 
order liquid-gas coexistence curve that culminates in a critical point p], and is 
extendable into the supercritical region as the Widom line [2], [3], jl] charac- 
terized by the locus of points with maximum correlation length £. Historically, 
such phase coexistence curves were modeled by the van der Waals (vdW) equa- 
tion augmented by the Maxwell "equal area" construction. This approach lies 
at the foundation of the modern thermodynamic picture characterizing coexist- 
ing phases through equal Gibbs free energies. However, the vdW-Maxwell theory 
suffers from several unresolved conceptual drawbacks [5], [6], and furthermore, an 
analytic prediction of the Widom line from any equation of state is yet unknown. 
Here we devise a novel construction to characterize liquid-gas phase transitions 
based on the continuity of £ between the phases, with the Riemannian geometric 
thermodynamic scalar curvature \R\ ~ £ 3 [7], which also allows, for the first time, 
a direct computation of the Widom line. The idea that the correlation lengths 
of the coexisting phases are equal, and its computational realization through R, 
provides a method for predicting the phase coexistence curve when used in con- 
junction with any equation of state, or experimental data. We illustrate this point 
here with the vdW equation, settling a century old problem of thermodynamic 
computation. We also determine the location of the Widom line for several fluids 
both with the vdW equation and with data from the NIST database [8]. Our 
results may be used to predict phase behaviour in a wide variety of systems, from 
boiling water to black holes, and promises to have significant impact on diverse 
areas of physical sciences and engineering. 

The key physical idea in our analysis originates from the microscopic picture 
of first-order liquid-gas phase transitions due to Widom [9]. In this framework, 
spontaneous density fluctuations cause the local density p (f) in a single phase 
fluid to deviate from the overall mean density po in some complex, time dependent 
manner. Mathematically, p(f) = po corresponds to an intricate contour surface 
that separates two sides with local mean densities p > po and p < p . A straight 
line through the fluid intersects this surface at points spaced an average distance 
£ apart, where £ is the correlation length characterizing the size of organized 
structures inside the fluid. £ is generally small in a disorganized system like an 
ideal gas, but diverges at the critical point for real fluids. When a single phase 
fluid is very near a first-order phase transition, a small amount of a second, 
minority phase will begin to form. A reference point in this single phase fluid 
typically has local density close to that of either of the two incipient coexisting 
phases. The typical density difference |Ap| across the contour surface p (f) = 
Po thus equals that of the two phases. Reversing the role of the majority and 
minority phases leaves this argument unchanged, with the same |Ap|. £ in the 
single phase plays a similar role in anticipating the properties in the two phases 
since £ is the thickness of the interface between the two [9]. This anticipated 
interface thickness must be the same approaching the phase transition with either 
of the two phases being the majority phase, and the correlation length £ should 
thus be the same in the two coexisting phases, the condition at the heart of our 
approach. 

For experimental predictability, we need an estimate of £, allowing a thermo- 
dynamic expression for the equality of the correlation lengths at the interface. 
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Figure 1: -R for the coexisting liquid and gas phases versus (T — T c ) /T c for 
normal Hydrogen calculated with the NIST fluid database, plotted on the loga- 
rithmic scale. At the indicated value, \R\ ~ v g , where v g is the molecular volume 
in the gas phase. Below this value of \R\, its interpretation as the correlation 
length loses significance. 

This can be realized using the Riemannian geometry of the equilibrium ther- 
modynamic state space of the system through the metric g a p = —j^q^§^f ( a > 
(3 = 1,2), where ks is Boltzmann's constant, and s, a 1 , and a 2 denote the en- 
tropy, energy, and particle number per unit volume, respectively [7J. The metric 
is based on Gaussian fluctuation theory whose breakdown takes place when the 
volume of the system is of the order of the Riemann scalar curvature R of the 
thermodynamic metric. This volume is expected to be the correlation volume 
£ 3 , leading to the desired connection [7J, \R\ ~ £ 3 . Experimental predictions 
for the coexistence curves of first-order phase transitions can thus be obtained 
from the equality of \R\ calculated in the two coexisting phases. We call this the 
i?-crossing method. At the critical point, R diverges. In the supercritical region 
beyond the critical point, the locus of the maximum of \R\, via \R\ ~ £ 3 , provides 
a theoretical prediction of the Widom line. 

As a direct test of our proposal, we calculate R for Hydrogen in both phases 
using the NIST fluid database [8], [10] and its program REFPROP. These pro- 
vide data based on phenomenological equations of state, with fit parameters 
determined by matching to experimental data for fluids. Results are shown in 
fig.(l), where agreement between the R's in the two phases is seen to be excellent 
near the critical point, better than 1% in the range 0.96 < T/T c < 1, where T is 
the temperature and T c its critical value. By contrast, at T/T c = 0.96, the molar 
densities of the coexisting liquid and gas phases differ from each other by a factor 
of ~ 3. 

Our i?-crossing method complements the canonical macroscopic rule for first- 
order phase transitions, namely the equality of the molar Gibb's free energies g 
of the coexisting phases pj. Applied to the vdW equation however, this macro- 
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Figure 2: R vs p r along an isotherm of Helium with t r = 0.86 in the lower half 
and v r vs p r along the same t r in the upper half (where the v r values have been 
multiplied by a factor of 3). The blue and green curves represent the stable 
branches, and the red curve is the unstable branch. We mark by "i?-crossing" 
the p r where the i?'s of the liquid and gas phases become equal (with c v = 1.5 
and 1.2 on the gas and liquid sides respectively). The line labeled "Maxwell" 
represents the corresponding p r from Maxwell's construction. 



scopic equality has unresolved conceptual problems. Finding states with equal g's 
involve contentious issues of integration along a reversible path through a thermo- 
dynamically unstable region in the Maxwell construction, or through the critical 
point in Kahl's approach [TT] . Such conceptual difficulties, which have been de- 
bated for over a century, are entirely bypassed in our construction. Further, our 
method naturally contains a measure of its limit of applicability, since for the 
construction to be effective, £ 3 should be large enough to encompass a number 
of atoms adequate for a thermodynamic approach to be reasonable. This limits 
us to a regime not too far from the critical point. We find that the i?-crossing 
method retains its viability down to volume regimes containing but about a single 
molecule. 

As a simple theoretical example, we apply the .R-crossing method to the uni- 
versal vdW equation in its reduced form, 



Pr 



8t T 



3v r — 1 



(1) 



where p r = P/P c , t r = T/T c 



v/v c and P and v are the pressure and molar 



volume, with the subscript c denoting their critical values. The critical quantities 
are known to be related to the vdW constants a and b by P c = T c = 27 ^° T 
and v c = 36. The Maxwell equal-area construction yields the limiting slope of the 
coexistence curve dp r /dt r = 4, independent of the fluid and its heat capacity. Our 
.R-crossing method inherits the same limiting slope here. This number is closely 
followed only by Helium and Hydrogen, for which this example is expected to 
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Figure 3: Phase coexistence and the Widom line for Helium (T c = 5.19K, P c = 
2.26 bars) on the left and Hydrogen (T c = 33.19K, P c = 13.30bars) on the 
right. The coexistence curve is calculated from vdW with the Maxwell equal-area 
construction and with i?-crossing, and compared with experimental data from 
NIST [8J. The Widom line is calculated by finding the locus of maximum values 
of \R\ both with vdW and from NIST data. We compare with the maximum 
values of c p from experimental data (NIST). The liquid and gas heat capacities 
c v i and indicated for vdW. In the supercritical region, we use c vg 



be maximally effective. R can be calculated here via the thermodynamic metric 
using standard formulae [7] and gives R = A ■ B, where 



A = ~- , 7V —2 > ( 2 ) 



3v 

3 c v (p r vf — 3v r + 2) 

and 



B = c v (p 2 r v 5 r - 9 Pr v? + 12 Pr v 3 r - 27v 2 r - Pr v 2 r + 27v r - 3) 

+18v r ( Pr vl + 1) (3) 

where c v is the dimensionless molecular specific heat at constant volume (assumed 
constant, though possibly different in the liquid and gas phases) and b plays no 
role in our subsequent analysis. 

For vdW isotherms with given t r < 1, substituting p r from eq.([T]) into eqs.(j2D 
and ([3]) results in two physical branches for R, corresponding to the liquid and 
gas phases (see color coded fig. (2)), with R diverging at the end points. The value 
of p r where the R values are equal (i.e they cross) is interpreted as the reduced 
saturation pressure corresponding to t r . For the cases we consider here, c v on the 
gas side is taken as 1.5, the ideal gas value. On the liquid side, we have chosen 
average values determined from NIST data [8], over the range of temperatures 
that we are interested in. Equivalently, for vdW isobars, the i?-crossing method 
can be used to predict the saturation temperature. 

In the supercritical region, isobaric \R\ exhibits a local maximum with respect 
to t r , whose locus is naturally interpreted as the Widom line, signifying a crossover 
for certain dynamical fluid properties from gas like on the low pressure side to 
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1.6 
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180.41 
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170.49 



Table 1: Saturation temperatures on the left and Widom line temperatures on 
the right (in Kelvins) for Neon (T c = 44.49 K, P c = 26.79 bars) on the upper 
part and Argon (T c = 150.69 K, P c = 48.63 bars) on the lower part. TJ, at (vdW) 
is the prediction of the saturation temperature from the i?-crossing method, us- 
ing the vdW equation, and is compared with experimental values from NIST. 
Corresponding values of \R\/v g are also shown to indicate the validity of our 
method. Widom line predictions from the i?-maximization method are obtained 
both from vdW with c v = 1.5 (T^(vdW)), and from NIST (T^f (NIST)). We 
have also shown the prediction of the Widom line obtained from maximising c p 
asT5. 



liquid like on the high pressure side [2], [3], [4]. We can calculate the Widom 
line as per its definition through \R\ ~ £ 3 ; free from the theoretical difficulty of 
characterizing it via the maximum of the specific heat Cp cLS IS conventional in the 
literature [3]. To find maxima for \R\ and c p , we search along isobars. 

A natural estimate for the validity of our analysis for vdW is offered by the 
dimensionless quantity \R\/v g , where v g is the coexistence molecular volume in 
the gas phase. \R\ > v g implies that £ 3 is greater than a molecular volume, and 
we are in a regime where our analysis is strictly valid. We find that for vdW, this 
restricts us to t r > 0.8 along the coexistence curve (a value indicated in fig. (3) 
and table (1)), and to p r < 10, in the supercritical region. 

Figure (3) summarizes our results for Helium and Hydrogen. Table (1) sup- 
plements these for Neon and Argon. From fig. (3), it can be seen that the R- 
crossing method, in conjunction with vdW, predicts excellent results within its 
range of applicability. Away from criticality, deviation from data is also due to 
the mean field nature of the vdW equation of state. Direct application of the 
i?-maximization method using NIST data in the supercritical regime shows strik- 
ing agreement with experimental c p maximum values even far from the scaling 
region. 

In conclusion, we have constructed a novel geometrical technique to charac- 
terize liquid-gas phase transitions from a microscopic perspective, through the 
thermodynamic scalar curvature R. When applied in conjunction with the vdW 
equation, this bypasses theoretical issues a century old with the Maxwell equal 
area construction and its variants. Our technique generalizes to any phenomeno- 
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logical equation of state, including those obtained as multi parameter fits to 
experimental data. This analysis further provides the first direct theoretical con- 
struction of the Widom line, without using any ad hoc thermodynamic response 
function. 

Our method unifies concepts in Riemannian geometry, thermodynamics, phase 
transitions, critical and supercritical phenomena. Although we have primarily 
applied our technique to liquid-gas phase transitions, the method should be uni- 
versally applicable to any first-order phase transition. This makes it of crucial 
significance to a diverse range of disciplines in physical, chemical and biological 
sciences, and engineering. It further generalises even to gravitational systems like 
anti de-Sitter black holes which also appear to exhibit liquid-gas like first-order 
phase transitions [T2] . 
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